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ON  TK£  PRINCIPX*E  OP  INVARIANT  IMBEDDING 
AND  DIFFUSE  REFLECTION  FROM  CYLINDRICAL  REGIONS 

By 

Richard  Bellman  and  Robert  Kalaba 

Introduction . 

The  Importance  of  varlouf3  principles  of  Invariance  In 

treating  problems  of  radiative  transfer  has  been  discussed  by 
12  4 

Ambarzumlan  ’  and  Chandrasekhar  ,  who  applied  these  principles 
to  a  variety  of  problems  connected  with  plane— para  1  lei  regions. 
However,  the  application  of  these  concepts  to  the  corresponding 
problems  for  cylindrical  and  spherical  regions  Is  not  Immediate. 
Cnandrasekhar^  has  observed:  "The  extension  of  the  approximate 
methods  of  solution  to  other  geometries  Is  straightforward;  but 
It  Is  not  equally  apparent  what  the  general  principles  are  which 
will  play  the  same  unique  role  as  the  principle  of  Invariance 
In  plane— paral lei  atmospheres." 

In  our  first  paper^,  we  introduced  the  principle  of  Invar¬ 
iant  Imbedding  as  an  extension  of  the  Ideas  of  Ambarzumlan  and 
Chandrasekhar,  and  utilized  this  principle  to  treat  the  case  of 
Inhomogeneous  plane— para  1  lei  regions.  In  this  paper,  we  propose 
to  show  that  this  new  concept  enables  us  to  treat  cylindrical,  and 
analogously  spherical,  geometries  In  a  unitary  fashion,  and  thus 
provides  an  answer  to  the  query  of  Chandrasekhar. 
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Anaiytlc  and  computational  questions  arising  from  this  treat¬ 
ment  will  be  discussed  In  greater  detail  In  a  subsequent  publication. 
2.  The  Physical  Model . 

Let  us  now  consider  the  problem  of  determining  the  Intensity 
of  diffusely  reflected  radiation  from  the  surface  of  a  circular 
cylinder.  We  assume  that  the  symmetry  properties  of  the  medium  are 
such  that  It  suffices  to  consider  Interactions  In  a  cross— sect  Ion 
perpendicular  to  the  axis,  and  that  the  absorption  and  scattering 
parameters  of  elemental  volumes  depend  only  upon  distance  from  the 
axis,  and,  finally,  that  the  scattering  Is  Isotropic. 

Consider  an  Incoming  flux  of  radiation  making  a  constant  angle 
0  with  the  radius  of  the  circle,  as  Indicated  In  the  figure  below, 
and  taken  to  have  Intensity  one  per  unit  length  normal  to  the  direc¬ 
tion  of  propagation.  This  Is  the  appropriate  definition  of  flux 
for  circular  regions  corresponding  to  the  usual  one  for  plane  regions. 


Reflected 

Radiation 


As  In  the  plane  case,  the  ratio  of  scattering  coefficient  to 
absorption  coefficient  ’  ?  denoted  by  A(x )/( 1— A(x )),  where  x  Is  the 
distance  from  the  center  of  the  circle. 

We  now  consider  the  Intensity  of  flux  of  reflected  radiation 
making  a  fixed  angle  ¥  with  the  extended  radius  of  a  cylinder  of 
radius  x,  due  to  an  Incoming  flux  at  angle  0,  as  above.  This 
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Intensity  Is  considered  to  be  a  function  r(v,u,x)  of  u  -  cos  0, 
V  ■  cos’f',  and  x.  In  this  way,  we  Imbed  the  particular  process 
within  a  continuous  family  of  processes. 

3-  Application  of  the  Principle  of  Invariant  Imbedding . 

In  order  to  obtain  a  functional  equation  for  r(v,u,x),  we 
consider  the  scattering  from  a  cylinder  of  radius  x  dx  as  the 
result  of  scattering  from  the  Infinitesimal  annulus,  (x,x-»-dx), 
and  the  cylinder  of  radius  x. 


Although  we  are  definitely  considering  multiple  scattering, 
our  functional  equations  are  considerably  simplified  by  the  fact 
that,  to  higher  order  terms,  only  one  scattering  can  occur  within 
the  annulus. 


Proceeding  as  In  Reference  3  which  In  turn  follows  the 
procedures  of  Ambarzumlan  and  Chandrasekhar,  we  obtain  the  fol¬ 
lowing  nonlinear  Integro-dlf f erentlal  equation 


together  with  the  boundary  condition,  r(v,u,o)  =•  0. 
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This  dlffere  from  the  corresponding  equation  for  the  Inho¬ 
mogeneous  paral lei— plane  case,  derived  In  our  first  paper^,by  the 
preeence  of  the  partial  derivatives  with  respect  to  u  and  v.  These 
terms  arise  from  the  fact  that  the  radiation  Incident  at  angle  9 
on  a  circle  of  radlua  x  +  dx  is  Incident,  to  first  order  terms,  at 
an  angle  x  (tan  9)  dx/x  on  a  circle  of  radius  x  — >  0. 

4  .  Reduction  Dimensionality . 

As  the  equation  stands.  It  Is  an  equation  for  a  function  of 
three  variables,  r(v,u,x),  two  space— like  variables,  u  and  v,  and 
one  time— like  variable,  x.  Let  us  now  show  that  this  equation 
may  be  replaced  by  two  equations  for  the  functions  P(v,x)  and 
Q(u,x)  defined  as 

1 

P(v,x)  -y^r(v,z,x)  dz/z, 

^  (1) 

Q(u,x)  -y^r(w,u,x)  dw. 
o 

The  equation  In  (5.1)  has  the  form 

L(r)  -  ?(P(v,x),  Q(v,x)),  (2) 

where  L  Is  a  linear  partial  differential  operator.  Solving  for  r 
In  terms  of  P,  using  characteristic  theory  or  otherwise,  we  obtain 
a  relation  of  the  form 


r 


A(') 


(?) 
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where  /\18  a  linear  operator.  Using  thla  relation,  we  form  P(v 
and  Q(u,x)  as  defined  above. 

1 

P(v,x)  •J'  A(P)  dz/z, 

0 

(4) 

1 

Q(u,x)  /\(p'  dw. 

o 

This  last  set  of  equations  constitute  a  system  of  two 
nonlinear  Integral  equations  for  the  functions  P(v,x)  and  Q(u,x 

Computation  and  analytic  aspects  of  these  equations  will 
be  presented  subsequently. 
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